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Abstract. Previous studies have found that calculations which consider long- 
range magnetic dipolar interactions truncated at a finite cut-off distance Rc predict 
spurious (unphysical) long-range ordered phases for Ising and Heisenberg systems 
on the pyrochlore lattice. In this paper we show that, similar to these two cases, 
calculations that use truncated dipolar interactions to model the Gd3Ga50i2 garnet 
antiferromagnet also predict unphysical phases with incommensurate ordering wave 
vector qQj.(j that is very sensitive to the dipolar cut-off distance Rc- 
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1. Introduction 

There are currently many highly-frustrated magnetic materials being experimentally 
studied where the magnetic species consist of a rare-earth 4f ion, such as Ho^+, Dy^+, 
Qj3+ 'j]-,3+^ which can have a large magnetic dipole moment. Because of the large 
moment, the long-range dipole-dipole interactions in these systems are an important 
part of the full spin Hamiltonian. The role of dipolar interactions in highly frustrated 
magnetic Ising systems has been systematically investigated for the three-dimensional 
pyrochlore lattice of corner- sharing tetrahedra [ll[2l|33. 

In highly-frustrated Heisenberg antiferromagnets of corner-sharing triangles or 
tetrahedra, any state with zero total magnetic moment on each elementary triangle 
or tetrahedron unit is a classical ground state [1]. There are an infinite number of 
such spin configurations and this is why these systems fail to develop conventional 
magnetic order at nonzero temperature [1]. In cases where dipolar interactions are 
somewhat weaker than nearest-neighbor exchange interactions, one might have naively 



assumed that the long-range dipolar interactions can be truncated at a finite cut-off 
distance, Rc, since the nearest-neighbor exchange energetically controls and enforces 
the nearest-neighbor correlations. Previous studies on the Ising spin ice pyrochlore 
systems [H [21 [3] and the Heisenberg pyrochlore antiferromagnet [5l [6] have found that 
this naive expectation is erroneous and that truncating the dipolar interactions at a 
finite cut-off distance Rc leads to spurious (unphysical) long-range ordered phases, and 
that it is crucial to consider dipolar interactions to infinite distance {Rc = oo). In 
this paper we report a third example of a highly frustrated spin system which is very 
sensitive to the dipolar cut-off. Specifically, we consider a Heisenberg model on the 
three-dimensional garnet lattice structure of corner-shared triangles. Our work extends 
the study of a dipolar Ising version of the model on a garnet lattice [7] and is relevant to 
the ultimate understanding of the nature of the incommensurate spin-spin correlations 
that develop in the Gd3Ga50i2 garnet (GGG) antiferromagnet below a temperature of 
500 mK [8l[9]. 

2. Model and method 

In order to investigate the problem of an adequate treatment of the dipolar interactions 
in GGG, we consider below a minimal model Hamiltonian H for it. To best expose 
the physics of a truncated dipolar lattice sum, we ignore the effects of the quantum 
nature of the Gd^"*" spins, lattice disorder [8], exchange interactions beyond nearest 
neighbors |H1 UHl [11] and possible single ion anisotropy [12] , all of which are potentially 
important for a thorough quantitative understanding of GGG. We describe the magnetic 
Gd^+ spins S as classical and isotropic n = 3 component (Heisenberg) vectors of the 
length ^ (S(S + 1)) (S = 7/2) coupled by frustrated antiferromagnetic nearest-neighbor 
exchange and long range dipolar interactions, of strength J = 0.107 K and D = 0.0457 
K, respectively |TOl |TT] : 

H= J 5] S, • S, + D 5: [S. • S, - 3 (S, ■ %)(S, ■ f,,)] . (1) 

<i,j> i>j «i 

In Eq. ([1]), i,j span the sites of the GGG lattice (see Fig. 2 in Ref. [10] for the GGG 
lattice structure) which are separated by vectors r^j = i^ij'f^ij of directions fjj; 
denotes pairs of nearest neighbors. The size of the GGG conventional cubic cell is 
a = 12.349 A|H]. 

We aim to identify the critical (or, soft) modes of model ([1]) for which a magnetic 
instability first develops as the temperature, T, is reduced. We consider the soft mode 
spectrum in the Gaussian (mean- field theory, or MFT) approximation [13]. Following 
Ref. [T3|, we apply MFT to calculate the neutron scattering intensity /(q). This allows 
for a convenient way of analyzing the physical infiuence of finite Rc on the magnetic 
correlations as well as for a direct comparison with experimental data [8]: 

m = [fmr x hm i/Nj:{st ■ s,^)e^''-- . (2) 




Figure 1. Each panel shows the upper branch A"P(q) of eigenvahies of model ^ 
as a function of wave vector q in the (hhl) plane at various cut-off distances Re- 
Rc — 3,4,5 (upper panels, left to right) and Rc = 100, 1000, oo (lower panels). The 
global maximum Amax defines the ordering wave vector qord (denoted by an arrow). 
Its location is a complex function of Rc- The isolines are drawn at 0.5%, 3%, 7%, 13%, 
of the overall dispersion downhill from the global maximum; the actual values of the 
dispersions as well as of Amax depend on Rc (q is measured in units of 27r/a). 

Here, angular brackets denote a thermal average, is the number of spins, represents 
the components of spin Sj at site i perpendicular to the scattering vector q, and /(|q|) 
is the magnetic form-factor of Gd^"*"!!!]. The MFT expression for /(q) [13] is obtained 
from the eigenvalues A"(q) and eigenvectors of the basis-diagonalization of the Fourier 
transform of exchange and dipolar interactions in Eq. ([T]) [13] : 

Here a = {1, . . . ,n ■ Nf, = 36} enumerates the eigenvalues and the vector F"(q) 
incorporates information on the eigenvectors and represents the role of the paramagnetic 



form factor of the GGG primitive unit cell, which contains Ni, = 12 ions. The ordering 
wave vector qord is given by locating in the first Brillouin zone the global maximum, 
Amax, of the maximum value (upper branch), A"P(q), among the 36 A"(q) eigenvalues. 
The mean field critical temperature, T^^^ = Xmax/n, is used to define a (positive) 
dimensionless temperature r = T/T^^"^ — 1 to serve as a natural energy scale. 

3. Results 

With the aim of studying the effect of the dipolar cut-off on the magnetic correlations of 
model ([1]), we first compute, for various Rc, A"P(q) for arbitrary q in the first Brillouin 
zone. Technically, we calculate the A"(q) modes on a finite 32^ q-space grid in the 
zone, and obtain their values at any q using a three dimensional cubic interpolation 
procedure. We verify the grid-independence of the results by considering denser grids 
and by cross-checking the interpolated results with exact calculations at judiciously 
chosen q values. 

We find that the dipolar term of Eq. ([T]) selects a unique ordering wave vector 
qord with corresponding mode A'^P(qord) out of the massively degenerate spectrum of 
soft modes of the nearest-neighbor model at any cut-off distance Rc > 1. The ordering 
wavevector qord was found to belong to the (hhl) planes of the first Brillouin zone for 
all Rc > 1. However, its location in those planes is very sensitive to Rc- We display this 
in Fig. [Hby showing the dependence of the A"P(q) modes on q in the (hhl) plane for 
Rc = 3,4,5, 100, 1000. We note their three important properties as reflected in Fig. [H 
First, for each value of Rc, A"P(q) is characterized by a relatively small overall dispersion 
Amax/Amin — 1 ~ 10% throughout the zone, where Amin is the global minimum of the 
branch. Even within the interval 1% below Amax (or at 13% of the overall dispersion, as 
delineated by the outermost isolines), A'^P(q) covers the major part of the (hhl) plane. 
Second, though the general topology of A'^P(q) within the Brillouin zone is preserved 
with varying dipolar cut-off distance Rc, its fine structure is manifestly sensitive to Rc- 
The ordering wave vector qord displays a non-monotonous dependence upon Rc as shown 
in Fig. [T]by its movement throughout the zone. Third, qord converges to a well-defined 
value, and so does the dispersion of A"P(q) away from qord, at large values of Rc, as is 
evident by comparing the Rc = 1000 panel with the limiting case of Rc = oo, recast by 
the Ewald method [13j. These are the central results of the paper. 

We now proceed to calculate the neutron scattering intensity /(q) of model ([1]) 
within the MET scheme, Eq. Unlike the spectra A"(q), /(q) can be directly 

compared to experiments, such as the one on a powder sample of GGG in zero external 
magnetic field [8]. We determine the powder intensity I{q) by numerically calculating 
the spherical average of /(q), which entails the application of a cubic interpolation 
procedure separately to the numerator and denominator of Eq. ([3]). To illustrate the 
influence of Rc on the spin-spin correlations of model ([H), we show in Fig. [2] the MET 
powder scattering profiles /(g) of model ([1]) at several Rc = 3, 10, 100, oo. 

The /(g) profiles show different degrees of dependence on the dipolar cut-off distance 




Figure 2. Theoretical powder neutron scattering profiles of model ([T]) at the 
dipolar cut-ofF distances Rc = 3, 10, 100, oo and the dimensionless temperatures 
T = 0.1, 0.01, 0.001, 0.0001. At r = 0.1 the profiles are not sensitive to the cut-off 
distance and reflect in fact the properties of the nearest-neighbor Heisenberg model 
without dipolar interactions. With the approach to the critical temperature the profiles 
acquire dispersions that are unique reflections of different ordering tendencies of model 
([T]) at different Rc- The profiles are uniformly off-set for clarity. 



i?c at different dimensionless MFT temperatures r (Fig. [2]). At temperatures sufficiently 
far from the mean-field critical regime (Fig. [2], r = 0.1), the magnetic correlations are 
not very sensitive to Rc- In fact, it has been found in Monte Carlo simulations [8], 
that paramagnetic liquid-like correlations of GGG can be well described by completely 
ignoring the dipolar term. Ai t 0.01 the profiles start to capture the effect of the 
dipolar interactions on the magnetic correlations. The effect turns out to depend on Rf.. 
Indeed, at r ~ 0.0013 the profiles clearly display a specific g— dependence sensitive to the 
chosen R^.. Starting from this regime, the influence of small Rc becomes uncontrollable, 
as seen by the formation of Bragg peaks at spurious ordering wave vectors (cf. Fig. [1]). 
Figure [2] shows that even the consideration of a rather large dipolar cut-off of Rc = 100 

t The powder MFT Bragg intensities grow as log |t| as opposed to the l/r growth of the q-dependent 
intensities. This explains the necessity of considering rather small r in order to theoretically approach 
the critical regime. 



does not allow one to reproduce magnetic correlations consistent with the physical 
Rc = oo limit. This, together with the incommensurability of the fundamental ordering 
wave vector of the physical Rc = oo case: qord = ^n/a (0.348 0.348 0.253) strongly 
warns against using a standard Monte Carlo method with periodic boundary conditions 
to tackle this problem. Moreover, we anticipate that a finite-size correction of 

the real space representation of the Ewald interactions would be sufficiently large for 
numerically accessible system sizes so as to prohibit a quantitative disentanglement of 
the role of perturbative terms PllTl|12] to model ([I]), presumably a necessary condition 
for obtaining a quantitative description of the experimental incommensurate magnetic 
correlations fHl in GGG. 



4. Conclusion 

To conclude, we have identified another example of a highly frustrated Heisenberg 
antiferromagnetic system, namely that on the garnet lattice, where the selection of 
the soft mode is sensitive to an ad-hoc cut-off distance Rc of the dipolar interactions. 
This adds to the cases of the Ising (spin ice) [H El E] and the dipolar Heisenberg 
antiferromagnets [5l [6], both on the pyrochlore lattice. Continued progress in 
understanding GGG may be possible provided the dipolar interactions are treated at 
their physical infinite cut-off limit Rc = oo. This will, however, require systematic 
investigations of the roles of long-range exchange [H] and single ion anisotropy [12] 
in this material. Finally, we note that the positions and relative intensity of the I{q) 
maxima differ dramatically from those found experimentally (Fig. 2a of Ref. [8]); an 
adjustment of r does not solve the discrepancy. This may indicate that a quantitative 
description of the low-temperature spin-spin correlations in GGG requires inclusion of 
exchange interactions beyond nearest neighbours flU[ ITT] and/or single ion anisotropy. 
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